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Abstract 

The question of the averaging of inhomogeneous spacetimes in cosmology is im- 
portant for the correct interpretation of cosmological data. In this paper we suggest 
a conceptually simpler approach to averaging in cosmology based on the averaging of 
scalars within unimodular gravity. As an illustration, we consider the example of an 
exact spherically symmetric dust model, and show that within this approach averaging 
introduces correlations (corrections) to the effective dynamical evolution equation in 
the form of a spatial curvature term. 
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1 Introduction 



The Universe is not isotropic or spatially homogeneous on local scales. The correct gov- 
erning equations on cosmological scales are obtained by averaging the gravitational field 
equations (FE). An averaging of inhomogeneous spacetimes in Einstein's general relativity 
(GR) can lead to dynamical behavior different from the spatially homogeneous and isotropic 
Friedmann-Lemaitre-Robertson- Walker (FLRW) model; in particular, the expansion rate 
may be significantly affected [T]. Consequently, a solution of the averaging problem is of 
considerable importance for the correct interpretation of cosmological data. 

The solution to this problem necessitates a method for covariantly (and gauge invariantly) 
averaging tensors on a background differential manifold. Unfortunately, this is a very difficult 
problem. In the Isaacson spacetime averaging scheme in macroscopic gravity (MG) bilocal 
averaging operators are utilized [2j. Choosing a compact region £ C M. in an (n- dimensional 
differentiable) manifold (Ai, g a p) with a volume n-form and a supporting point x G S to 
which the average value will be prescribed, the average value of a geometric object, Pp(x), 
over a region E (with volume Vs) at x G S, is defined in terms of the bilocal extension of 
the object Pp(x), Pp(x, x') = W™,(x, x')p^,(x')Wp (x' , x), by means of the bilocal averaging 
operator Wp,(x,x'). The averaging scheme is covariant and linear by construction, and the 
averaged object has the same tensorial character as pg. In any manifold with a volume 
n-form there always exist locally volume-preserving divergence- free operators j2], in which 
the bilocal operator (x', x) takes the simplest possible form: Wp'(x',x) = SfSp [3]. 

The definition of an average consequently takes on a particularly simple form when 
written in a volume-preserving (system of) coordinates (VPC). Indeed, if the manifold is 
a pseudo-Riemannian spacetime, the spacetime average of a tensor field Pg(x), x G £, at a 
supporting point (t, x a ) G £ in VPC is thus 

( P «(t,x a )) e = ^rj Pp(t + t\x a + x a ')dt'd 3 x' . (l.i) 
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In the MG covariant approach to the averaging problem the Einstein FE (EFE) on 
cosmological scales with a continuous distribution of cosmological matter are modified by 
appropriate gravitational correlation (correction) terms j2]. The averaged FE can always 
be written in the form of the FE for the macroscopic metric tensor when the correlation 
terms are moved to the right-hand side of the averaged field equations, and consequently 
can be regarded as a geometric modification to the averaged (macroscopic) matter energy- 
momentum tensor [2j. In [3j it was found that by solving the MG equations the averaged 
EFE for a spatially homogeneous, isotropic macroscopic spacetime geometry has the form of 
the EFE of GR for an FLRW geometry with an additional spatial curvature term (i.e., the 
correlation tensor is of the form of a spatial curvature term) (see also [5]). Unfortunately, the 
spacetime averaging scheme in MG is very difficult to apply and is fraught with complications 
[B]. In this paper we want to suggest an alternative approach to averaging, exploiting the 
preferred nature of VPC and based on the averaging of scalars [H |8]. 



2 Unimodular gravity 



The fundamental variables in the action for unimodular gravity and the Einstein-Hilbert 
action for GR are different [9J [TOJ E] . In unimodular gravity, there is an additional restriction 
on the metric, not present in GR: the determinant of 9 „„ e q nal S one. A s a consequence of 
det g^ v = 1, unimodular gravity is only invariant under volume-preserving diffeomorphisms. □ 
Thus, unimodular gravity presents a natural theory in which to do averaging. 

Varying the action in unimodular gravity leads to the FE relating the traceless Ricci 
tensor, Rjf = — ^Rg^, to the corresponding traceless energy-momentum tensor T£ v [9]. 
It should be noted that the energy-momentum conservation law V M T M1/ = does not follow 
from this equation of motion, but has to be imposed separately. Assuming energy-momentum 
1 Coordinate invariance can always be reinstated into the theory. 
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conservation, it then follows that 



R + T = -A (2.2) 

where A is a constant (and = 1 and c = 1). Using the contracted Bianchi identity, 
VuG^ u = 0, we then obtain 



where A is given in terms of A and the vacuum energy density (part of the energy-momentum 
tensor) p vac = A vac . Hence the cosmological constant A naturally appears in terms of a 
constant of integration in unimodular gravity. 

Therefore, the theory acquires a new integrability condition [9]. Both the trace-free FE 
and the matter conservation equations are assumed; the integrability condition follows from 
these equations. Hence, we obtain the differential relations which are functionally equivalent 
to the full EFE (where the spacetime volume density yjg is not a dynamical variable), where 
the cosmological constant is thus given in terms of an arbitrary integration constant A and 
is not given explicitly by the vacuum energy A vac . 

This is an old proposal essentially initiated by Einstein [12] and more recently it has 
been developed under the name of unimodular gravity [T0| [T3| ITT] . It has been suggested 
that unimodular gravity can be used to eliminate problems caused by the nature of the 
cosmological constant as well as to resolve the discrepancies between theory and observation, 
while not introducing any exotic terms such as quintessence or dark energy into the analysis 
of the EFE [HJ [JUJ [TT] . Indeed, although unimodular gravity does not give a unique value 
for the effective cosmological constant, it has the potential to solve the huge discrepancy 
between theory and observation. With a suitable high-energy cut-off, the vacuum energy 
density is estimated by Weinberg [9] to be of the order p vac ~ 2 x 10 71 GeV 4 , whereas the 
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effective value of the cosmological constant as determined by astronomical observations is of 
the order p Q ^ s ~ 10 -47 GeV 4 . However, there is no longer a cosmological constant problem. 
For example, for a perfect fluid the matter source term is the manifestly trace-free stress 
tensor (p + p) (u a Ub + (1/4) g a b)'i hence, matter enters the FE only in terms of the inertial 
mass density (p + p), which vanishes in the case of a cosmological constant (e.g., see |15j). 

Unimodular gravity has also been utilized in the study of the quantization of GR [TTl[T3] . 
The Hamiltonian of a generally covariant theory is zero, so in a sense there is no evolution, but 
since unimodular gravity is not generally covariant, the classical problem of time is avoided 
[TT] . In addition, in unimodular gravity quantum gravitational factor ordering ambiguities 
are alleviated [T3] . 

3 Averaging Proposal 

We wish to exploit the structure of unimodular gravity to suggest an alternative approach to 
averaging in cosmology. Within unimodular gravity we need to average the trace-free part 
of the FE and the trace of the FE separately. 

(1) Average trace-free part of the FE: Here the resulting correlation tensor must consequently 
be trace-free. If the form of the resulting equations are of the algebraic form of a 'perfect 
fluid', as in the cosmological application (with a large scale FLRW geometry), then the 
correlation tensor must be of the form of an effective energy momentum tensor T e ^j for 
which the trace T e ^ = —p + 3p = 0, corresponding to a radiation fluid |16j . Note that if 
the matter is dust, then T tot = T dust + T e ^ = -p d - p r + 3p r = -\p d + 2p r ] + [p r + 3p r ], 
which could be (trivially) reinterpreted as a renormalized dust term (with energy density 
[pd + 2p r ]) and a term corresponding to a constant spatial curvature (with [p r + 3p r ]) 011]. 

(2) Average the trace of the FE: In this case we only need to work with the (generalized) 
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3.1 



Friedmann eqn. 

The problem of averaging is then effectively reduced to considering the average of a single 
scalar eqn. (see [?]). 



4 Example: Lemaitre-Tolman-Bondi model 



The exact spherically symmetric dust Lemaitre-Tolman-Bondi (LTB) model [18] , which can 
be regarded as an exact inhomogeneous generalization of the FLRW solution, can be rewrit- 
ten in VPC (t,x,u,4>) [5]. Taking A = A(t,x), the line-element becomes 

alu 2 



als 2 



U 2 \ , 2 U , , dx 2 a2 
— \ d t>-2— A dtdx+ — + A 2 



+ (1 - u l )d^ 



(4.3) 



which has det(g) = 1 as desired, where U(t,x) is defined as 

u <m = - 2 t*+A A " - < 4 - 4 > 

The constraints on the original LTB metric ensuring a dust solution are given in [5]. For 
general functions A and U, the Ricci scalar of the metric (14. 3 j) is given by 



i2 = (1 - 5^A 4 + 3t/ 2 ^ - 2U x AA t - 2A XX A 5 + 3A 2 + 6UA x A t 

-2AUU X A X + A 2 UU XX + A 2 U tx + A 2 U 2 ) (4.5) 



The spatially flat (E = 0) FLRW model in VPC is given by the metric (14. 3 j) with 

2 Note that the sum of the averaged energy-momentum tensor and the correlation tensor is covariantly 
conserved; the question of whether the averaged energy-momentum is separately conserved with respect 
to the averaged geometry is determined by averaging the energy-momentum conservation equation (if it is 
not, then there is an effective interation between the averaged energy-momentum tensor and the correlation 
tensor [T7]). 
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A(t,x) = A = (3x) 



1/3 



U(t,x) 



(4.6) 



t-t B 7 



where, strictly speaking, the degenerate form for U(t, x) does not follow directly from eqn. 
(1P|) (however, eqn. fl4~5j) is valid for (146]) and R ~ (t - t B y 2 ). Defining A = rS(t), 
the Ricci scalar of the FLRW metric with positive curvature constant k = +1 is given by 
R = Q[SS tt + S t 2 + k]S~ 2 . For the zero-curvature Einstein de-Sitter metric, S ~ and 
setting t B = 0, we get the approximate expression: 



consistent with the expression given in [5] (with E = \k = cr 2 ). 
4.0.1 A Perturbative Solution 

let us assume that ts(r) is zero, which implies that the bang time is uniform and we are 
consequently restricting our choice of LTB models to those with no decaying modes. We 
shall also consider solutions of the LTB metric in VPC as perturbations about the spatially 
flat FLRW model given by f)4.6p . In this respect our approximate solution will be an expan- 
sion with respect to E and we require the Einstein tensor to have the form of dust (after 
truncation of terms of O(E^) or higher). We begin by making the formal expansion for A 
in the form: 



where a±, «2, a, b, c and d are constants. We can use eqns. (14. 8 h and (I4.4p to obtain 
U(t,x). Calculating the Einstein tensor and requiring it have the form of dust (up to order 
Eq) allows us to determine the constants in our perturbative solution (we obtain: a = 1/3, 
b = 0, c = 5/6 and d= -1 0). 

The expression that results from substituting U in terms of A using equation (I4.4p and 



R=^t- 2 [l + htl + 0(t)}, 



(4.7) 



A(t, x) = A + axX^Eo + a 2 x c t d E 2 , 



(4.8) 
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the expression (14.81) for A (with the given powers of x and t in our particular perturbative 
solution) leads to the expression for the Ricci scalar R (keeping only terms up to O(Eq)): 

R= — - AE a lX - 2 / 3 - El (Ih ai x- 1 IH- 1 + 2a 2 2 ^x- 2 ^) . (4.9) 
3t z 

Defining r 3 = 3x t -2 , we obtain 

R = + aE oai r~ 2 t-^ + E 2 (ba^r 4 / 3 + cajr-H^ 3 ) (4.10) 

where 

a = -4 x 3 2/3 , b = -15 x 3 1/6 , ce-2x 3 1/3 . (4.11) 

Finally, we obtain the averaged version of the Ricci scalar equation by integrating eqn. 
(14. 10p over the radial variable r, where ro is the (radial) averaging length scale: 



R = JL (i + aE oai t 2 / 3 + E 2 [b ai t 2 / 3 + colt 2 ' 3 ] ) (4.12) 

Ob 

(where the 'barred' constants are the appropriately r -renormalized constants). We see that 
all of the correction terms (correlations) introduced by averaging the Ricci scalar equation 
are of the form of a spatial curvature term (j4.7p . which is consistent with the results of [5]. [] 



5 Discussion 



Recent observations are usually interpreted as implying that the Universe is very nearly 

flat, currently accelerating and indicating the existence of dark matter and dark energy 

3 We note that for this perturbative solution, we obtain higher order correction terms of the form ~ t~ 2 , 
which can be interpreted as a renormalization of Aq in the exact dust solution. We also note that, in 
principle, for the second order terms (0(£q)) to be formally comparable, a\ 2 ~ o^o'- 
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[T9~] . A cosmological constant is a candidate for the dark energy. Averaging can have a very 
significant dynamical effect on the evolution of the Universe; the correction terms change the 
interpretation of observations so that they need to be accounted for carefully to determine 
if the models may be consistent with an accelerating Universe. Indeed, it has been argued 
that a more conservative approach to explain the acceleration of the Universe without the 
introduction of exotic fields might be to utilize a backreaction effect due to inhomogeneities 
of the Universe. 

In this paper we have argued that a rigorous approach to cosmological averaging (and 
necessary for studying cosmological data) is perhaps most naturally studied within the con- 
text of unimodular gravity. In the simple example studied here, we found that all correction 
terms introduce correlations of the form of a spatial curvature term [5]. 

As another simple illustration, in the special case (C e = 0, e« = 0) of the exact solu- 
tion representing a two-scale Buchert average of the EFE for an inhomogeneous universe 
approximating the observed Universe [20], we have that S(t) = at 2//3 [l + flt] 1 / 3 , where 
a = a (3# ) 2/3 (l - /,o) 1/3 (2 + f v0 )- 2/3 and /3 ee 3/ u0 #o(1 - + f v0 )'\ and the 

physical constant f V Q is related to the scales representing the voids and the bubble walls 
surrounding them (within which clusters of galaxies are located). In this example, the Ricci 
scalar is again of the form of eqn. (14.71) . 

In future work we intend to consider this averaging scheme in more general cosmological 
contexts. In particular, we wish to study approximate solutions within linear perturbation 
theory. A first step will be to develop perturbation theory within unimodular gravity |21j. 
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